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SUMMARY 


Analytically-derived  plane  strain  stresses  and  strains  in  an  annulus 
pressurised  sufficiently  to  cause  plastic  flow  are  given.  Unloading  giving  rise 
to  reyielding  around  the  bore  is  then  examined,  along  with  the  effect  of 
reaming.  Cold-working  is  then  included  in  the  analysis  in  terms  of 
interference  between  mandrel  and  hole.  Finally,  comparisons  of  stress  and 
strain  predictions  are  made  with  those  from  a  finite  element  analysis.  ■R.Vv3'Thcx\  tO 
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Superscript 


2  In  (c/a)  +  1  -  (c/b)2 
Lame  constant  (equation  4) 
change  (applied  to  radius) 
strain 

dimensionless  interference,  I /a 
Poisson’s  ratio 
reyield  radius 
stress 


at  bore 

at  elasto-plastic  radius 

mandrel 

yield 
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circumferential 
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1.  V  INTRODUCTION 

u 

The  fatigue  process  in  aircraft  structures  is  very  often  associated  with  the 
presence  of  stress  concentrations  of  one  kind  or  another,  and  fastener  holes 
represent  by  far  the  most  numerous  such  feature  in  any  mainly  metallic  structure. 
For  some  time  now  substantial  attention  has  been  given  to  methods  for  improving  the 
fatigue  behaviour  of  structures  containing  holes,  one  of  the  most  successful  of  these 
methods  being  the  cold-expansion  process’^"  In  common  with  other  mechanical 
means,  hole  cold-expansion  achieves  its  quite  dramatic  effect  (life  improvement 
factors  typically  of  the  order  of  two  to  five)  by  modifying  the  local  stress  field*'  to 
become  benign  so  far  as  the  fatigue  process  is  concerned. 

In  the  cold-expansion  process  the  hole  is  expanded  by  the  passage  of  an 
oversize  or  interference-fit  tapered  mandrel  through  the  hole,  a  sleeve  or  protective 
•  bush  being  used  to  prevent  the  possibility  of  scoring  or  galling^  As  the  mandrel 
r^^entenTthe  hole,  the  surrounding  material  is  first  expanded  (typically  of  the  order  of  3 
to  5%  of  diameter)  causing  local  yield  to  occur.  As  the  mandrel  leaves  the  hole  the 
surrounding  elastically  deformed  material  forces  a  reduction  in  hole  diameter  from 
the  fully  expanded  size,  and  this  reduction  can  be  such  as  to  cause  plastic  flow  in  the 
reverse  or  negative  sense.  Of  course  the  extent  (i.e.  radius)  of  this  reyielded  zone  is 
substantially  less  than  that  generated  by  the  mandrel  in  enlarging  the  hole.  The  final 
stage  in  cold-working  is  usually  the  reaming  of  the  hole  to  provide  a  precise  diameter 
to  receive  a  transition-  or  interference-fit  fastener.  It  can  be  seen  that  the  process 
results  in  a  region  of  material  surrounding  the  hole  which  is  compressively  stressed, 
and  this  is  highly  beneficial  in  delaying  the  onset  of  fatigue  crack  initiation  and  in 
reducing  the  rate  of  propagation  of  developing  fatigue  cracks.  It  is  not  until  the 
growing  crack  tip  reaches  that  region  corresponding  to  a  tensile  circumferential 
stress  field  that  the  beneficial  effect  of  the  cold-working  process  is  lost  -  this  occurs 
typically  at  about  2  to  3  hole  radii  -  and  by  that  time  the  sought  substantial  increase 
in  total  fatigue  life  has  already  been  achieved. 

The  analysis  of  this  report  relates  to  an  axisym metric  thick-walled  annulus 
under  conditions  of  plane  strain,  pressurised  and  depressurised  to  simulate  the  cold- 
expansion  process.  The  elasto-plastic  pressurised  cylinder  has  been  much  studied  in 
the  past  4_1C>  and  this  report  makes  no  further  contribution  to  this  field  per  se.  It 
does,  however,  focus  the  application  of  this  work  on  pressurised  cylinders  to  hole 
cold-expansion  in  annuli.  Use  is  made  of  Hencky  deformation  theory,  incompressible 
plastic  flow  and  the  von  Mises  yield  criterion,  all  of  which  have  previously  been 
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applied  to  the  elasto-plastic  cylinder  problem  .  A  recent  paper  on  cold-expansion 
also  used  these  same  assumptions.  To  the  author’s  knowledge,  however,  no  single 
source  provides  formulations  of  all  of  the  quantities  of  interest  relating  to  cold- 
expansion:  it  is  for  this  reason  that  the  present  report  has  been  prepared. 

For  the  sake  of  convenience  and  completeness  here,  the  derivations  of  all 
stresses,  strains,  elasto-plastic  radii,  etc  are  given  for  each  of  the  three  main  cold- 
expansion  stages  as  follows: 

(a)  mandrel  fully  in  -  hole  fully  expanded,  plastic  flow  around  hole 

(b)  mandrel  removed  -  effect  of  elastic  recovery,  reyielding  at  hole 

(c)  effect  of  reaming  -  modification  to  stress  and  strain  fields 

Then,  because  cold-expansion  is  specified  in  terms  of  interference  between  mandrel 
and  hole,  rather  than  applied  pressure,  the  relations  between  these  for  the  various 
conditions  are  derived  to  permit  recasting  of  the  stress  and  strain  equations  in  terms 
of  interference. 

Finally,  a  comparison  is  made  between  stress  and  strain  predictions  from  the 
present  report  with  those  from  a  plane  strain  finite  element  study  using  incremental 
plasticity  theory.  The  two  sets  of  predictions  show  good  agreement  for  the 
quantities  of  interest. 

2.  GENERAL  RELATIONSHIPS 

In  the  axisymmetric  problem  to  be  discussed,  the  radial,  circumferential  and 
axial  directions  are  also  the  principal  directions.  All  formulations  are  given  in  terms 
of  these  principal  directions. 

2.1  Elastic  Region 

For  plane  strain  conditions  (.e^  =  0)  the  following  general  relationships 

hold5: 

dur 

er  =  Car  -  v  (o9  +  V]/E  =  dT 

ee  =  loe  -  V  (oz  +  or)l/ E  =  -JT1 

6  =  Ca  -  v  <o„  +  a  )]/E  =  0 

z  z  r  9 


(1) 
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From  the  last  of  these 


o,  =  v(cj  +  o„) 

z  re 


(2) 


Rearranging  (1)  to  obtain  stresses  in  terms  of  strains  gives 


*r  =  y  Ier  +  l^£e1 


°e°r[‘e+  *  «r3 


(3) 


°Z  =y[er+  ee3 


where  y  is  given  by 


y  =  E 


1-v 


(l+v)(l-2v) 


(4) 


2.2  Plastic  Region 

A  substantial  and  extending  literature  exists  on  the  intricacies  of  how  to  best 
represent  plastic  flow.  All  theories,  of  necessity,  include  approximations  which 
better  or  worse  suit  a  given  application.  _  this  report,  for  the  sake  of  achievability 
in  conjunction  with  practical  usefulness,  deformation  plasticity  theory  is  used.  The 
following  assumptions  are  made; 

(a)  Incompressibility  -  plastic  flow  occurs  with  no  change  in  volume.  This 
corresponds  to  Poisson’s  ratio  becoming  equal  to  one  half. 

(b)  Plane  strain  -  as  before,  6Z  =  0>  and  from  the  Hencky  equations6  it 
follows  that 

°z  =  %  (or  +  °0>  (5) 

and 


€r  +  6e  =  0 


(6) 
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These  results  also  follow  directly  from  (1)  and  (2)  with  v  = 

(c)  Yield  criterion  -  in  terms  of  principal  stresses,  the  von  Mises  criterion 
is  given  by 

<or-  o0)2  +  ^e-az>2+  <az-  or>2  =  2a02 
Substituting  for  azfrom  (5)  gives 


(d)  An  elastic/perfectly-plastic  stress-strain  relationship  is  assumed  in 
both  positive  and  negative  senses.  Work  hardening  does  not  occur  and 
the  Bauschinger  effect  is  not  accounted  for. 

Because  of  these  simplifications,  incompatibilities  do  occur  -  in  particular, 
because  of  the  need  for  step  changes  in  the  values  of  v  across  elasto-plastic 
boundaries.  For  practical  purposes,  however,  these  inconsistencies  are  usually  either 
small  in  magnitude  or  affect  quantities  of  secondary  interest  only. 

3.  STRESSES 

In  this  section  the  principal  stress  distributions  are  given  for  the  annulus  first 
under  internal  pressure  loading  and  then  following  removal  of  the  pressure. 

3.1  Loaded  Annulus 

The  application  of  pressure  to  the  bore  of  the  annulus  gives  rise  to  three 
possible  outcomes  of  interest: 

(a)  elastic  expansion  throughout 

(b)  expansion  such  that  the  material  at  the  bore  is  about  to  yield  - 
incipient  yielding 

(c)  expansion  such  that  a  plastic  zone  extends  from  the  bore  to  an  elasto- 
plastic  boundary,  the  surrounding  material  remaining  elastic. 

Each  of  these  is  now  examined  in  turn. 


V 
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3.1.1  Elastic  Conditions  Tlirougfaout 

The  equilibrium  equation  for  an  annulus  under  axisymmetric  loading  is  given 


by" 


°r  -  n 
dr  +  r 


(8) 


the  solution  of  which  is5,  for  internal  pressure  loading  p, 


°r  = 


-  P 


(r>  ~  (P 


=  P 


r>  + 


1  -  <r> 


(9) 


and  from  (2), 


2  up 


By  way  of  example,  these  stress  distributions  for  the  case  of  b/a  =  5 
and  v  =  1/3  are  shown  in  Fig.  2.  Note  that  a z  is  constant  across  the  section. 


3.1.2  Incipient  Plastic  Flow 

Increase  in  pressure  in  the  bore  eventually  leads  to  the  onset  of  plastic  flow 
there.  Substitution  of  (9)  into  (7)  provides  the  pressure  at  which  plastic  deformation 
begins: 


1  -  (a/b)2 

pa=  °°  VC3  +  (l-2v)2(a/b)“] 

where  the  subscript  denotes  incipient  plastic  flow  at  r  =  a. 

For  practical  values  of  a/b  in  the  cold-expansion  situation  <no  larger  than 
one-half),  the  second  term  in  the  denominator  above  is  quite  negligible,  so  that,  for 
practical  purposes 


P 


=  -77  [1  -  <r>  ] 
a  v3  b 


(10) 
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3.1.3  Stresses  in  Plastic  Region 

Furthe*  increase  in  pressure  gives  rise  to  an  enlarging  region  of  plastic  flow 
around  the  hole,  surrounded  by  a  restraining  elastic  region  beyond.  The  elastic 
region  provides  a  constraint  on  the  extent  of  plastic  flow  which  is  crucial  to  the 
effectiveness  of  the  cold-expansion  process  in  increasing  fatigue  life. 

In  the  plastic  region,  as  already  noted  at  (7),  the  yield  criterion  becomes 


a 


e 


(11) 


Substitution  of  (11)  into  the  equilibrium  equation  (8),  integrating,  and  noting 
that  at  r  =  c  the  radial  stress  there  is  given  by  (10)  with  radius  a  replaced 
by  c,  leads  to 

r  r 2 

O  --  -T  [2  In  7  -  1  +  <£>J 

r  v3  c  b 


From  (11) 


"6  -73  [2  ln  c  +  1  +(b)] 


(12) 


and  from  (5) 


In  7  +(£)] 
c  b 


Equations  (12)  give  the  stress  distributions  in  the  plastic  region  (a  s  r  s  c> 
of  a  pressurised  annulus  surrounded  by  an  elastic  region  between  radii  c  and  b. 

The  pressure  (applied  at  r  =  a  )  to  cause  yielding  to  radius  c  is  found  from 
the  first  of  (12)  as 


V  <-Vc  ■  t![2  ,n  ‘7’  +  1  -  ‘f’1  <l3> 

O'o 

=  73  ° 


where 


a  =  2  In  f  +  1  -  (£j 
a  b 


(14) 
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3.1.4  Stresses  in  Elastic  Region 

The  stresses  in  the  elastic  region  are  found  from  (9)  by  substituting  for  p  the 
expression  (10)  and  by  changing  all  a  to  c.  Thus 


0  c2  c2 

*r  -  V3  ^  +  <f)] 


=  7J  C  <r > 


+  (b)] 


(15) 


2  vo 

~V3 


0  c2 
<r> 


Equations  (15)  give  the  stress  distributions  in  the  elastic  region  ( c  s  r  s  b ) 
of  a  pressurised  annulus  in  which  the  region  a  s  r  s  c  has  undergone  plastic  flow. 
Note  that,  as  before,  o z  is  constant  in  this  region,  and  for  b  large,  ■*  0. 

The  elasto-plastic  stress  distributions  across  the  section  for  the  same  case  as 
before,  and  with  b/c  =  2,  are  shown  in  Fig.  3.  This  shows  the  discontinuity 
in  at  r  =  c  caused  by  the  step-change  in  v  at  the  elasto-plastic  boundary  -  it 
decreases  with  increasing  b  and  disappears  when  b  becomes  infinite.  More 

n 

sophisticated  and  more  complicated  analyses  introduce  second-order  terms  to 
eliminate  this.  In  the  context  of  this  report,  third-direction  parameters  are  given  for 
completeness  only. 


3.2  Unloaded  Annulus 

Upon  removal  of  the  pressure  in  the  bore,  recovery  of  the  material  takes 
place  under  the  action  of  the  outer  elastically  stressed  material.  The  elastic 
recovery  is  quantified  by  the  superposition  of  the  elastic  equations  (9)  with  an 
appropriate  sign  change.  The  extent  of  the  unloading  is  governed  by  the  requirement 
that  the  residual  radial  stress  at  the  bore  of  the  hole  must  be  zero.  In  general,  this 
gives  rise  to  another  (reyielded)  plastic  zone  in  the  neighbourhood  of  the  bore,  of 
much  smaller  radius  than  that  occurring  under  load.  Alternatively,  since  the 
pressure  applied  to  the  hole  in  Section  3.1  is  given  by  (13),  this  is  precisely  the  value 
required  to  be  substituted  for  p  in  (9)  to  establish  the  residual  stress  fields  by 
subtraction  of  (9)  from  (12)  and  (15). 
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Three  outcomes  are  possible: 

(a)  elastic  unloading  across  the  full  width  of  the  annulus. 

(b)  as  above,  but  with  incipient  reyielding  at  the  bore. 

(c)  as  above,  but  with  a  finite  reyielded  region  surrounding  the  bore. 

It  is  assumed  in  the  following  that  loading  of  the  annulus  has  given  rise  to  a 
fully  elasto-plastic  situation  as  detailed  in  Sections  3.1.3  and  3.1.4  (Fig.  3).  For  the 
elastic  and  incipient  plastic  loaded  situations  of  Sections  3.1.1  and  3.1.2  only  elastic 
unloading  is  possible,  leaving  no  residual  stresses. 

3.2.1  Elastic  Recovery  -  No  Reyielding 

When  the  recovery  from  the  loaded  configuration  is  wholly  elastic,  no  further 
plastic  flow  is  induced  at  the  bore.  In  this  case  the  residual  stresses  are  found  as 
outlined  above. 


For  the  zone  as  r  s  c 


a  =  "77  <C2  In  -r  -  1  +  +  a 

r  V3  c  b 


<7>  -  <{> 

_r _ d 

a  2 

1  -  (b> 


2  >  + 

-  (C2  In  7  +  1  +  (f)  ]  -  «  — - 77“  > 


J9  V3 


*  -  f> 


(16) 


z  V3 


r  ,c 
c  b 


<f> 


1  "  (b> 


For  the  zone  c  s  rsb 


°r  -  73  +  <b>  ]  +  “  -7-7 if 

1  ~  'b; 


°e  =  V3  <r>  +  (P  1  "  * 


(7>  +  <b> 

-<t> 


(17) 


-  9  - 


°z  = 


2vo0  r2 

7T  <  (f  > 


<r> 


1  -  <£> 


where,  as  noted  at  (14), 

«  =  2  In  7  +  1  -  (Jo 
a  b 

Equations  (16)  and  (17)  give  the  residual  stresses  in  an  annulus  in  which 
loading  has  caused  yield  to  radius  c . 


3.2.2  Elastic  Recovery  -  Incipient  Reyield  at  Bore 

In  this  case,  since  =  0  at  the  bore,  the  yield  criterion  (7)  becomes 


CT 


e 


(18) 


Substituting  for  o  from  (18)  in  the  second  of  (16),  and  putting  r  =  a  there 

t 7 


gives 


"  *  <C2  m  ?  +  1  +  (b)  3  "  a 


1  * 

^-<t> 


3> 


from  which 


2  In  f  +  1  -  <r>  =  2  Cl  -  (f-)] 
a  b  b 


(19) 


Equation  (19)  specifies  the  condition  under  which  reyield  will  initiate  at  the  bore  on 
unloading  in  terms  of  the  elasto-plastic  radius  c . 

Rearranging  (19)  as 


<r>- 


<f>  -  2 

O 

2  In  7  -  1 

a 


(20) 


and  solving  for  b/a  provides  the  locus  of  incipient  reyielding  shown  in  Fig.  4. 
For  b/a  less  than  2.22  reyielding  can  never  occur,  and  reyielding  will  never  occur 
when  c/a  is  less  than  1.65.  The  first  condition  is  a  geometric  requirement,  the 
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second  a  loading  requirement:  there  has  to  be  sufficient  surrounding  material 
available  to  impose  the  necessary  recovery  on  unloading,  and  the  extent  of  the  load- 
induced  plastic  zone  has  to  be  such  that  reyield  can  occur  on  unloading. 
For  b/a  greater  than  22.2,  reyielding  always  occurs  from  a  fully  plastic  expansion, 
i.e.  from  the  condition  c/a  =  b/a. 


3.2.3  Unloading  giving  rise  to  Reyielded  Zone 

In  this  case  unloading  gives  rise  to  a  reyielded  zone  around  the  bore  of 
radius  p,  as  indicated  diagramatically  in  Fig.  1(b). 


For  the  zone  a  s  r  s  p 


Stresses  in  the  reyielded  region  are  readily  found.  Substitution  of  (11)  into 
the  equilibrium  equation  (8),  integrating,  and  noting  that  =  0  at  r  =  a  provides 
the  following 


a 


r 


a 


e 


2°o  r 

-7TC1  +  ln7] 


(21) 


«2  =  -Sn+21nI3 


For  the  zone  parse 

Stresses  in  this  region  may  be  written  down  directly  from  (16)  after  first 
making  the  substitution  (19)  for  a  and  then  replacing  a  by  p: 

°r  '  S  <C2  >"  c  '  1  +  <b’]  +  2  [<f>  -  <b>]> 


o  =  -r  <C 2  In  7  +  1  +  (7)]  -  2  [<*)  +  (£>]> 
e  V3  c  b  r  b 


V3 


jr 

c 


c2 


(22) 
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In  this  region  reaches  a  minimum  when  r/p  =  V2,  at  which  point 

o  -a.  .  This  is  also  the  region  in  which  o  changes  sign:  that  occurs  when  r 
r  B  “ 

satisfies 


In 


(23) 


For  the  zone  c  i  r  s  b 

Stresses  in  this  region  may  also  be  written  down  directly.  From  (17),  after 
substituting  (19)  for  a  and  replacing  a  by  p : 


°r  =  73  +  <b)]  +  2  [(r>  '  (b)3> 


■  73  {C(r)  +  (b)]  -  2  C<r>  +  (b)]> 


(24) 


°°  C2  n2 

O  -  Zv  C(r)  -  2  (f)3 
z  V3  b  b 

Equations  (21)  to  (24)  contain  the  unknown  reyield  radius  p .  Radius  p  is  not 
independent,  but  depends  upon  the  extent  of  prior  yielding,  c,  and  the  geometry  of 
the  annulus.  Its  magnitude  is,  however,  readily  found. 


3.2.4  Relationship  between  Yield  Radii  c  and  p 

The  results  of  the  previous  Section  provide  a  ready  means  of  establishing  the 
relationship  between  c  and  p.  By  equating  like  stresses  (or  or  oQ)  from  (21) 
and  (22)  at  r  =  p  there  results 

2  In  +  !  -  (f)  =  2  [!-(£)]  (25) 


This  the  reyield  equivalent  of  the  incipient  reyield  case  of  Section  3.2.2  and 
reduces  to  (19)  when  p  *  a.  Equation  (25)  provides  the  sought  relationship  between 
p  and  c.  As  before,  (25)  may  be  rearranged  as 


(*■>  - 
P 


(->  -  2 
_ 


In  —  —  -  ! 


p  P 


(26) 


-  12  - 


The  plot  of  (26)  is  shown  in  Fig.  5,  which  is  an  extension  of  Fig.  4  for  the 
incipient  reyield  case.  The  reyield  regime  now  has  the  additional  information  shown 
in  terms  of  loci  of  the  parameters  p/a  and  c/p.  It  is  seen  that,  for  a  given 
geometry  b/a,  as  the  yield  radius  c  (or  c/a)  increases  -  corresponding  to 
increasing  load  at  the  bore  in  the  first  place  -  the  generated  reyield 
radius  p  (or  p/a  )  also  increases,  as  expected.  Also,  increasing  c/a  gives  rise  to 
increasing  c/p,  the  ratio  of  the  reyield  radii. 

For  b  large,  from  (26) 


p  =  V(ac)/e^ 

For  the  loaded  example  considered  earlier  and  shown  in  Fig.  3,  the  unloaded 
equivalent  may  now  be  evaluated.  By  solving  (26)  for  b/a  =  5  and  c/a  *  2.5  (or  by 
scaling  from  Fig.  5)  the  corresponding  p/a  =  1.19,  giving  c/p  =  2.10.  Equations 
(21)  to  (24)  may  now  be  evaluated,  the  plots  being  shown  in  Fig.  6.  Once  again  the 
inconsistencies  in  z-direction  stresses  across  the  elasto-plastic  boundaries  are 
apparent.  Notice,  however,  the  significant  extent  of  the  circumferential  stress  field 
below  zero  (from  (26)  this  extends  to  r/a  *  2.01)  and  the  size  of  the  reyielded 
region  (p/a  =  1.19)  :  both  aspects  play  a  major  part  in  delaying  and  then  slowing 
fatigue  crack  growth  from  the  hole. 

3.3  Effect  of  Reaming 

For  engineering  reasons  it  is  common  to  ream  holes  after  cold-expansion:  the 
final  radius  of  the  hole  thereby  becomes  slightly  larger  than  that  immediately  after 
cold-working.  The  effect  of  this  on  the  stress  fields  may  be  considered  as  follows: 
the  reaming  process  clearly  exerts  no  influence  on  the  loaded  situation  -  all  loaded 
stresses  remain  as  before.  But  consider  now  that  the  bore  of  the  loaded  annulus  is 
enlarged  so  that  on  unloading,  recovery  occurs  on  the  basis  of  this  larger  radius:  this 
is  exactly  the  situation  existing  after  reaming.  The  effect  is  found  simply  by 
substituting  a'  (=  a  +  sa)  fora  in  the  unloading  components  of  the  equations*. 


*  Care  is  required  in  tracking  radii  here.  In  evaluating  effects,  the  reamed 

radius  should  be  compared  with  the  post-cold-expanded,  pre-reamed  value  (see 
Section  5).  The  effect  of  reaming  must  be  calculated  on  this  small  difference  in 
radii,  not  on  the  larger  difference  between  final  ream  radius  and  initial,  unloaded 
radius. 
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As  anticipated,  reaming  marginally  relaxes  the  desirable  effects  of  the  residual 
stress  distributions  effected  by  cold-working. 

3.3.1  Reaming  where  no  Reyield  has  occurred 

This  case  is  that  of  all-elastic  unloading.  From  Section  3.2.1  the  required 
solutions  are  found  by  replacing  radius  a  in  (16)  and  (17)  with  a 1  in  the  elastic 
unloading  components,  i.e.  those  with  the  a  multiplier. 

For  the  zone  a 1  &  r  s  C 


'r  ‘  V3 


ra!: 


i  I  2 


r  c2  (r>  "  (b> 

oL  =  —  <C2  In  £  -  1  +  <§)3  +  a'  — 1 - 


1  ~  (b> 


o'  =  ~r  <C2  In  7  +  1  +  (7)3  -  a' 
6  V3  C  D 


A  1  2  A  1  2 

(7)  +  (f> 

r _ d 


I  2  ^ 


l-(b> 


(27) 


o'-  ~  <C2  In  7  +  (7)]  -  2ua' 
z  V3  CD 


a  1 2 

<t> 


1  -  (t) 


TT> 


For  the  zone  c  s  r  s  b 


=  73  <C-(r>  +  +  «’  T  \a:-2> 

1  'b' 


<*>*-  & 
..  _ E _ e. 


o'e  =73  <[(r>  +  (b)] 


A  1  2  A  1  2 

(7>  -  (F>  , 

A  1  2  ' 

1  -  <b) 


(28) 


I  1  2 


CVOc 


<r> 


0  c2  b 

a'  =  — r-  {(£)  -  a*  - TT> 

z  V3  b  .  fa_7 

1  V 


where,  by  analogy  with  (14), 


"  2  1nf'+  1  -  (b> 


(29) 
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Equations  (27)  and  (28)  give  the  residual  stress  after  reaming  to  radius  a1  in 
an  annulus  in  which  loading  caused  yield  to  radius  c . 


3.3.2  Reaming  -  Incipient  Reyield  at  Radius  r  =  a  ' 

Following  the  procedure  of  Section  3.2.2  there  results  an  identical 
relationship  to  (19)  where  a  there  becomes  a1 : 

2  In  J,  +  1  -  (J>  =  2  Cl  -  (|-)3  =  ct '  (30) 

Thus,  the  information  contained  in  Fig.  4  applies  also  to  the  reamed  hole 
provided  the  change  in  radius  is  noted. 


3.3.3  Reaming  where  Reyield  has  occurred 

In  this  case,  in  addition  to  the  change  in  radius  from  a  to  a',  reaming  gives 
rise  to  a  change  (increase)  in  the  reyield  radius  p  to  p 1 :  the  (reyielded)  stresses  at 
both  original  and  new  bores  remaining  the  same,  Fig.  7.  The  new  stresses  are  found 
directly  from  the  results  of  Section  3.2.3. 

For  the  zone  a 1  i  r  s  p 1 


2cto  r 
V3  ln  a' 

2o0  r 

7TC1+  ln7'] 


°z  = 


uo  r 

73  C1  +  2  1n  a,] 


For  the  zone  p  *  s  r  s  c 


(31) 


o 


a 


a 


;  -  ^<t2  ,n7  - 1  *  'f'1  ♦  2  'f’31 


.  I  2 


'  =  <12  ln  7  +  1  +  (J)]  -  2  [f)2+  (£>]> 

Q  V3  c  d  r  d 


.  1  2 


0  r  c*  o 1 2 

1  =  -=  H2  In  7  +  (7)]  -  2v  <£)  > 
Z  CD  D 


(32) 
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For  the  zone  c  i  r  s  b 


0  c2  c2  o'2  o'2 

o'  =  -75-  <[-(7)  +  <r>3  +  2  [(f)  -  <^>  ]> 
r  %/3  r  d  r  b 


0  r  2  r2  n *  2  n12 

O*  =  77  {[ (7)  +  <7>]  -  2  [(%  +  (^>  ]> 
e  V3  r  b  r  b 


(33) 


0  C2  n  *  2 

o'  =  -r  2i/  {(7)  -  2  (f)  > 
z  V3  b  b 


3.3.4  Effect  of  Reaming  on  Reyield  Radius 

Following  the  procedure  of  Section  3.2.4,  there  results  a  similar  relationship 

to  (24): 


2  In  J77  +  1  -  <§>  «  2  Cl  -  (£)']  (34) 

As  before,  the  information  of  Fig.  5  holds  in  this  case  provided  that 
a  and  p  are  replaced  by  a '  and  p 1  . 

The  effect  of  reaming  on  the  change  in  reyield  radius  p  is  readily  found.  By 
equating  the  expressions  for  ( c/b ) 2  from  (25)  and  (34)  for  the  pre-and  post-ream 
situations  it  is  found  that 


In 


For  b  large. 


a.' 

a 


and  for  7  =  1 
« 


6a  . 
+  and 
a 


*  =  1 

p 


bp  _  1,  6a 

p  a 


(35) 


(36) 


(37) 


i.e.  a  small  increase  in  radius  caused  by  reaming  increases  the  reyielded  radius  by 
about  half  as  much. 

This  effect  can  be  seen  in  the  previous  example  reworked  here  with  an 
assumed  (and  enormous)  10%  ream.  The  circumferential  stress  distributions  before 
and  after  reaming  are  shown  in  Fig.  7(a),  the  detail  in  the  reyield  region  being  shown 
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enlarged  in  Fig.  7(b).  It  is  clear  that  the  effect  of  this  degree  of  reaming  on  the 
residual  hoop  stress  is  small:  at  more  realistic,  much  smaller,  degrees  of  reaming  the 
changes  will  be  quite  inconsequential.  Any  possible  effect  of  the  slightly  less 
beneficial  stress  distribution  on  the  fatigue  process  from  this  point  of  view  would 
have  to  be  judged  negligible. 


4.  STRAINS 

In  this  Section  principal  strains  for  the  various  loading  conditions  are  given. 
The  format  follows  that  for  stresses.  In  both  elastic  and  plastic  regimes  in  two- 
dimensional  radial  flow  strains  and  displacements  are  very  simply  related: 


e 


r 


(38) 


and 


e 


9 


(39) 


so  that  radial  displacement  is  found  directly  from  circumferential  strain  as: 


Ur  =  (40) 

For  this  reason,  strains  alone  are  explicitly  formulated  in  the  following. 

4.1  Loaded  Annulus 

As  before,  the  application  of  pressure  to  the  bore  gives  rise  to  the  three 
outcomes  of  interest: 

(a)  elastic  expansion  throughout 

(b)  as  above,  but  with  incipient  yield  at  the  bore 

(c)  elasto-plastic  expansion 

Strains  in  all  regions  are  readily  found  from  (1)  or  (6)  in  conjunction  with  the 
appropriate  stress  substitutions. 

4.1.1  Elastic  Conditions  Throughout 

The  substitution  of  (9)  into  (1)  provide  the  required  equations  of  elastic  strain 
in  a  pressurised  annulus.  They  are: 
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and 


R 

E 


(1  ±  v) 
3  2 

1  "  (b> 


{(-)  - 
r 


(1 


-  2u)(J)> 


R 

E 


1  + 


1  "  <b> 


Up 


+  (1  -  2v >  <*■)> 


(41) 


4.1.2  Incipient  Plasticity 

At  this  point  (10)  provides  the  value  of  pressure  p  to  just  initiate  yielding  at 
the  bore.  Substitution  of  (10)  into  (41)  gives 


f>  « -  <>  -  *»>«**> 


and 


(42) 


1  +  v 
V3 


<*>  <  i 


+  (1  -  2v)(J)> 


4.1.3  Strains  in  Plastic  Region 

Under  conditions  of  plastic  plow,  as  already  noted  at  (6),  the  sum  of  the 
principal  strains  is  zero.  Substitution  of  (38)  and  (39)  into  (6)  provides  the 
differential  equation 


(43) 


the  solution  of  which  is 


u 


r 


(44) 


where  C  is  a  constant. 

From  (39)  and  (44) 


(45) 


Equating  (42)  and  (45)  across  the  elasto-plastic  boundary  (found  by  putting  r  =  c  and 
byreplacing  a  by  c  in  (42)),  C  may  be  evaluated,  leading  to 


1  +  v  °  c2  c2 

er  --vPr  (f)C1  +  <!  -  2u>(f>3 


and 


(46) 


1  +  u  °°  C2  C2 

<=  =  on  +  a  -  2v)(^>] 

e  V3  E  r  b 

since  in  the  plastic  region. 

4.1.4  Strains  in  Elastic  Region 

Strains  here  are  written  down  directly  from  (41),  substituting  for  p  the  value 
given  by  (10): 


1  4 

V3 


up  - 


(1 


-  2v)(|)3 


(47) 


e 


e 


1  +  V 

V3 


0  c2  c 

r[(r}  +  (1  ~  2v){b 


)] 


A  small  inconsistency  will  be  noted  here  arising  from  the  fact  that  the 
assumed  incompressibility  in  the  plastic  region  whereby  e  =  -  e  does  not  carry 

i  y 

over  into  the  elastic  region.  Consequently,  radial  strain  at  the  elasto-plastic 
boundary  calculated  from  (46)  does  not  equate  with  that  from  (47)  at  the  same 
position.  The  difference  is  proportional  to  twice  the  second  term  in  the  square 
brackets  which  is,  in  general,  relatively  small.  For  v  approaching  one  half,  or  for 
b  large,  this  term  becomes  insignificant  and  the  discontinuity  vanishes. 

Figure  8  shows  the  circumferential  and  radial  strain  distributions  for  the 
example  already  given,  highlighting  the  small  discontinuities  in  radial  strain  at  the 
elasto-plastic  boundary  r/a  =  2.5. 
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4.2  Unloaded  Annulus 

As  before,  three  outcomes  of  interest  are  possible: 

(a)  elastic  unloading  throughout 

(b)  as  above,  but  with  incipient  reyielding  at  the  bore 

(c)  as  above,  but  with  reyielding 

Also  as  before  it  is  assumed  here  that  an  elasto-plastic  situation  exists  prior 
to  unloading. 

4.2.1  Elastic  Recovery  -  No  Reyielding 

Subtraction  from  (46)  and  (47)  of  elastic  strains  given  by  (41)  with  p 
replaced  by  (13)  provide  the  sought  solutions. 


For  the  zone  a  s  r  s  c 


1  +  v 

er  "  ~  v/3 


1  +  v  °° 


c2 
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c2 
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For  the  zone  c  s  r  s  b 
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c2 
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a 
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$ 
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b  i 

-  a 
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$ 
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b  i 

•  u 
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1  -  <i> 

where,  from  (14), 


a  =  2  In  7  +  1  -  (7) 
a  d 


(48) 


(49) 
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4.2.2  Elastic  Recovery  -  Incipient  Reyield  at  Bore 

The  situation  here  is  similar  to  that  above,  except  that  now  some 
simplifications  occur  because  of  (19),  whereby 

a  =  2  In  -  +  1  -  (£)  =  2  tl  -  (f)3 
a  o  o 

For  the  zone  a  a  r  s  c 

1  j.  „  a°  r1  r 2  a2  a2 

er  =  -  <<r)C1  +  <1  "  2i/)<^>]  -  2  C(^)  -  (1  -  2v)<J)3> 


(50) 


1  +  u  0  C2  C2 

{(S.)[i  +  (i  _  2v)(r-)] 
V3  E  r  b 


2  [(-)  +  (1 
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2v)(^)]> 


For  the  zone  c  s  r  s  b 


er  =  " 
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V3  E 


—  (C(^)  -  (1  -  2v)(J)]  -  2  Up 


(1  -  2u)(~)]> 


(51) 
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<c<r)  +  <i  -  2y)<r)3  -  2  C(r>  +  <1  - 


2i/)(J>]> 


4.2.3  Unloading  giving  rise  to  Reyielded  Zone 

As  before  the  reyield  zone  radius  is  denoted  by  p. 


For  the  zone  a  i  r  s  p 

In  this  reyielded  region  (44)  and  (45)  again  hold.  By  equating  (45)  and  (50) 
across  the  (reyield)  elasto-plastic  boundary  (found  by  putting  r  =  p  and  by 
replacing  a  by  p  in  (50))  the  constant  in  (45)  may  be  evaluated.  The  results  are: 

er  =  -  -jr  «7>tl  +  (1  -  2v)(|)3  -  2  p)  Cl  +  (1  -  2v)(^)3> 
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and 


(52) 
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For  the  zone  p  &  r  S  C 

Strains  in  this  zone  follow  directly  from  (50)  where  a  is  replaced 

by  p. 

1  +  v  °°  c2  c2  D2  D2 

e  =  —  7-  +  <1  -  2v ) (t) 3  -  2  C (■£•)  -  (1  -  2v)(^)3> 

r  V3  E  r  b  r  b 


(53) 
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V3  E 


For  the  zone  c  s  r  s  b 

As  above,  strains  here  are  given  by  (51),  after  substituting  p  for  a : 


1  +  l,  0  C2  T2  n2  n2 

^3—  J-  <C(^)  -  (1  -  2u)(^)]  -  2  [(^)  -  (1  -  2iO(J)3> 
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V3 
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As  before,  Fig.  8  shows  the  unloaded  strain  distributions  across  the  annulus 
for  the  previous  example.  Now  two  discontinuities  in  radial  strain  are  apparent,  one 
at  each  of  the  elasto-plastic  boundaries. 


4.3  Effect  of  Reaming 

As  already  discussed  under  Section  3.3,  the  effect  of  reaming  (enlarging  the 
bore  radius,  after  cold  expansion,  from  a  to  a 1 )  is  readily  found  by  replacing  a 
by  a'  in  the  unloading  terms  in  the  strain  equations  of  Section  4.2.  The  changes  in 
strain  (pre-  and  post-reaming)  are  compared  below  for  the  three  basic  cases 
considered  earlier. 


4.3.1  Elastic  Recovery  -  No  Reyielding 

The  strain  equations  for  the  pre-reamed  case  are  given  by  (48)  and  (49),  and 
for  the  reamed  case  by  the  same  equations  where  a'  replaces  a.  By  taking  the 
ratio  of  the  two  cases  the  effect  of  reaming  may  be  seen.  Unfortunately  the  forms 
of  (48)  and  (49)  do  not  lend  themselves  to  any  simplification  when  expressed  as 
ratios,  except  for  the  case  when  b  is  large.  In  that  case 
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independent  of  radius.  The  fractional  change  in  strain  is  the  same  in  both  radial  and 
circumferential  directions. 


4.3.2  Elastic  Recovery  -  Incipient  Reyield  at  Bore 

For  this  case  (50)  and  51)  apply,  and  for  b  large, 


at  any  radius. 


(Va' 

(Va 


Va' 

<Va 


(56) 


4.3.3  Unloading  giving  rise  to  Reyielded  Zone 

For  this  case  (52),  (53)  and  (54)  apply.  For  b  large,  and  at  any  radius 
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The  effect  of  the  ream  considered  in  the  example  of  Section  3.3.4  on  strains 
using  strain  ratios  from  (52),  (53)  and  (54)  is  shown  in  Fig.  9.  For  comparison,  the 
result  for  the  ordinate  using  (57)  is  0.82,  so  that  finite  annulus  geometry  must  be 
allowed  for. 


5 .  COLD  WORKING  -  EFFECT  OF  INTERFERENCE 

The  formulations  presented  in  Sections  3  and  4  for  stresses  and  strains 
contain  the  unknowns  of  pressure  (in  the  case  of  all-elastic  situations)  and  the  elasto- 
plastic  radii  c  and  p.  In  practice  the  magnitudes  of  these  parameters  arise  from 
the  interference  between  mandrel  and  hole  diameters.  In  this  Section  the 
relationships  between  interference  and  pressure  and  interference  and  elasto-plastic 
radii  are  developed. 
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In  the  cold-working  process,  as  the  mandrel  passes  through  the  hole,  the 
interference  originally  specified  is  accommodated  by 

(a)  a  radial  expansion  of  the  bore  of  the  hole  and 

(b)  a  radial  contraction  of  the  surface  of  the  mandrel. 

Thus 


I 


<u)  "  (um>a 

3  m  3 


(58) 


or 


<u)  ~ 

a  m  a 


A 


(59) 


where  A  is  the  dimensionless  (radial  or  diametral)  interference  and  the  subscript  m 
refers  to  the  mandrel.  Since  circumferential  strain  is  related  to  radial  displacement 
by  e0  =  u/r, 


(e 


e}a 


and 


A 


(Va  "  (<«)a 


(60) 


at  radius  r  =  a.  Since  the  circumferential  strain  for  the  annulus  has  already  been 
found,  it  remains  only  to  determine  that  for  the  mandrel  in  order  to  be  able  to  solve 
(60). 

Mandrels  are  used  in  practice  for  some  considerable  time  before  wear  takes 
its  inevitable  toll  -  typically  this  corresponds  to  cold-expanding  some  thousands  of 
holes.  During  this  time  they  remain  elastic,  and  their  analysis  in  terms  of  stresses 
and  strains  is  simply  that  of  an  externally  pressurised  solid  elastic  cylinder. 


5.1  Loaded  Mandrel 

The  mandrel  is  represented  as  a  pressurised  elastic  solid  cylinder  in  plane 
strain.  The  principal  stresses  are  given  by5 
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and 


(61) 


°z  =  '  2um  p 


all  stresses  being  independent  of  radius.  Circumferential  strain  is  found  directly 
from  (1)  by  substitution  of  (61)  as 


e 


em 


£■  (1  +  v  HI  -  2v  ) 
E_  m  m 


(62) 


Thus,  the  circumferential  strain  is  constant  at  all  points  throughout  the  mandrel, 
including  the  surface  at  radius  a,  the  displacement  there  being  given  by 

(u  )  =  -  £  (1  +  v  HI  -  2v  )a 

m  a  E„  m  tn 

m 

5.2  Loaded  Annulus 

The  three  basic  cases  already  considered  are  re-examined  below  in  terms  of 
interference. 


5.2.1  Elastic  Conditions  throughout 

Equations  (41)  for  circumferential  strain  represents  this  case.  Substituting 

(41)  and  (62)  into  (60)  gives,  with  Em  expressed  in  terms  of  the  modulus 

ratio  g  =  E/E  ,  at  r  =  a 
m 


£  _ 
E  ” 


A  (1  -  (|>) 

(1  +  v)Cl  +  (1  -  2 v)  (f)]  +  (1  +  vm)(l  -  2v  )  g  (1  -  (£>) 

b  m  m  b 


(63) 


Equation  (63)  defines  the  interface  pressure  in  terms  of  the  geometric  and 
material  properties  only.  Used  in  conjunction  with  (9)  and  (41),  stresses  and  strains 
in  the  annulus  may  be  evaluated  directly. 

In  situations  involving  interference-fitting  of  a  fastener  in  a  hole  it  is  normal 
that  the  plate  material  remains  elastic  -  interference  levels  are  usually  not  such  as 
to  cause  plastic  flow.  The  stresses  and  strains  in  such  situations  are  fully  described 
by  (9)  and  (41)  in  conjunction  wuh  (63)  when  both  fastener  and  annulus  are  in  a  state 
of  plane  strain. 


5.2.2  Incipient  Plasticity 


The  pressure  to  initiate  yield  at  the  bore  of  the  hole  has  been  found  (10)  as 

°°  a2 

P  -  v3  "  '  ‘b’1 

Substituting  this  value  for  p  in  (63)  gives 

A  —  =  {(1  +  v)Cl  +  (1  -  2v)(f)3  +  (1  +  v  )(1  -  2v  )  g  Cl  -  (£>]>  (64) 

o  q  \/3  d  m  m  b 

Equation  (64)  specifies,  for  given  geometry  and  material,  the  level  of 
interference  which  will  just  produce  yield  at  the  bore  of  the  hole.  Thus,  (64) 
determines  whether  plastic  flow  will,  or  will  not,  occur  for  a  given  level  of 
interference. 

For  the  case  of  common  Poisson's  ratios  and  moduli,  v  =  and 

g  =  1 ,  (64)  simplifies  to 


Thus,  the  onset  of  yield  is  independent  of  geometry  and  is  a  function  only  of  Poisson’s 
ratio  for  this  particular  case  only. 

5.2.3  Plastic  Region  -  Relationship  between  A  and  c 

In  the  plastic  region  equations  (46)  apply.  Substituting  (46)  and  (62)  into  (59) 
gives,  at  r  =  a , 

A  =  * ~  (7>tl  +  (1  -  2v)(£)3  +  £  (  1  +  V  >(1  -  2i'  )  (65) 

v3  t  a  D  m  m 

m 

The  pressure  p  applied  at  radius  r  =  a  to  cause  yielding  to  radius  c  has 
already  been  found  (13)  as 

°°  c  c2 

p  .  ^  [2  in  f  +  1  -  (f» 

Substituting  this  value  for  p  in  (65),  putting  =  E/g  and  rearranging  gives 
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A  —  =  77  {'1  +  v)<7)Cl  +  (1  -  2v)<7>]  + 
o0  V3  a  D 

+  (1  +  vm)(l  -  2um)  g  £2  In  J  +  1  -  (f)3>  (66) 

Equation  (66)  provides  the  sought  solution  between  interference  A 
and  elasto-plastic  boundary  radius  c:  unfortunately  it  is  clearly  not  practical  to 
express  c  in  terms  of  A.  The  effect  of  interference  on  the  elasto-plastic 
radius  c  as  a  function  of  geometry  b/a  is  given  in  Figs.  10  for  the  case 
v  =  v  =  1/3  :  in  Fig.  10(a)  the  modulus  ratio  g  =  1/3  and  in  Fig.  10(b),  g  =  1. 
Comparison  of  the  two  sets  of  graphs  shows  that  the  effect  of  g  is  not  marked,  and 
the  same  may  be  said  of  the  effect  of  geometry  -  the  plastic  radius  resulting  from  a 
given  interference  is  not  greatly  affected  by  b/a.  The  common  value  of 
AE/o0  occurring  at  incipience  (c/a  =  1)  for  the  case  g  =  1,  noted  in  Section  5.2.2, 
may  be  seen  in  Fig.  10(b).  Once  interference  is  such  as  to  locate  the  point  of 
interest  below  the  cluster  of  curves,  full  plastic  flow  has  occurred  throughout  the 
annulus.  This  is  a  situation  not  usually  found  in  cold-working. 

For  the  example  considered  earlier,  where  C/b  =  ^  ,  b/a  =  5  and 
c/a  =  2.5,  with  v  =  =  1/3  and  g  =  1/3  (typical  of  a  steel  mandrel  in  an 

aluminium  alloy  annulus)  Fig.  10(a)  (or  equation  (66))  gives 

AE/a0  =  5.43 

With  o0/E  =  0.005  for  a  typical  aluminium  alloy.  A,  expressed  as  a 
percentage,  becomes 

A  =  2.7%. 

Thus,  in  this  example,  an  interference  of  2.7%  gives  rise  to  an  elasto-plastic 
boundary  at  one-half  of  the  outside  radius.  For  a  more  representative  interference 
of  perhaps  3%  to  4%,  the  radius  of  that  boundary  will  be  a  little  larger. 

For  b  large,  (66)  becomes 

A  —  =  -7  {(1  +  v)(~)  +  (  1  +  um)(l  -  2v  >  g  [In  (7)  +  1]> 
a0  v3  a  m  m  ’  a 


(67) 
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The  first  term  is  dominating,  so  that  on  a  log-log  plot  a  virtual  straight  line 
results:  this  is  seen  in  Figs.  10. 

Finally,  it  should  be  noted  that  the  relationship  between  the  loaded  and 
unloaded  yield  radii  c  and  p  has  already  been  established  (25)  as 

2  In  ^7  +  1  -  <£>  *  2  tl  -  <f)3  (25) 

pz  D  D 

Thus,  for  a  given  interference,  c  may  be  found  from  (66)  or  Figs.  10,  and 
then  p  from  (25)  or  Fig.  5.  All  stresses  and  strains  in  plane  strain  cold-worked 
annuli  may  then  be  determined  from  the  appropriate  equations  in  Sections  3  and  4. 

6.  COMPARISON  WITH  FINITE  ELEMENT  ANALYSIS 

Finite  element  analyses  have  been  carried  out  under  two-dimensional  plane 
stress  and  plane  strain  conditions13  and  compared  with  results  from  a  three- 
dimensional  analysis14  in  which  hole  diameter  was  equal  to  the  thickness  of  the 
annulus,  inner  and  outer  radii  and  expansions  being  identical.  That  comparison 
showed  that  the  situation  in  the  interior  of  the  three-dimensional  annulus  was  very 
well  described  by  the  two-dimensional  results.  It  is  of  interest,  therefore,  to 
ascertain  how  well  results  from  the  formulations  in  this  report  agree  with  the  finite 
element  plane  strain  predictions. 

The  finite  element  anlyses  were  carried  out  using  the  incremental  plasticity 
routines  included  in  ARL’s  PAFEC  computer  codes.  Apart  from  the  differing 
theoretical  bases  (incremental  versus  deformation  plasticity  theory)  the  finite 
element  analyses  used  the  bilinear  work-hardening  stress-strain  relationship 
indicated  below,  whereas  the  present  report  is  based  on  an  elastic /perfectly-plastic 
relationship.  In  each  case  the  radius  ratio  b/a  =  10  and  a  4%  interference  between 
the  steel  mandrel  and  the  hole  in  the  aluminium  alloy  annulus  has  been  used. 


Property 

Annulus 

Mandrel 

Modulus  of  elasticity  (MPa) 

69,000 

209,000 

Poisson’s  ratio 

0.33 

0.30 

Yield  stress  (MPa) 

480 

Work  Hardening  slope  (MPa) 

1,200 

(finite  element  analysis  only) 
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6.1  Plane  Strain  Analysis  (this  report) 

In  any  given  case  it  will  usually  be  known  whether  plastic  flow  is  to  be 
expected.  If  in  doubt,  the  determining  equation  (64)  may  be  used.  In  the  present 
case  the  value  of  the  right  side  of  (64)  is  0.6%,  so  that  an  applied  interference  A  of 
4%  will  certainly  cause  yield.  Reference  to  Fig.  10(a)  or  the  solving  of  (66)  by  trial 
and  error  for  c/a  gives 

7  =  2.645 

a 

Reference  to  Fig.  5  indicates  that  reyielding  is  to  be  expected  and  (25)  may 
be  used  to  solve  precisely  (also  by  trial  and  error)  for  p .  Here 

£  =  1.255 

a 

The  several  radius  ratios  required  in  solving  the  stress  and  strain  equations 
are  then  noted.  These  are: 


J  =  0.265 

£  =  1.255 

a 

-  =  2.108 
p 

£  =  0.126 

and  the  required  solutions  are  then  calculated  from  the  following  equations  for  the 
cases  indicated. 

For  the  loaded  case 

Radial  Zone  Stress  Strain 


a  s  r  s  c  (12)  (46) 


c  s 
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For  the  unloaded  case 


Radial  Zone  Stress  Strain 


a  s  r  s  p  (21)  (52) 

p  s  r  s  c  (22)  (53) 

c  i  r  s  b  (24)  (54) 


6.2  Comparison  of  2-D  Solutions 

Loaded  and  unloaded  stress  and  strain  predictions  from  the  finite  element 
solutions  and  those  calculated  as  indicated  above  are  compared  in  Figs.  11,  12  and 
13. 

It  is  seen  that: 

(a)  circumferential  and  radial  stresses  show  good  agreement  for  both 
loaded  and  unloaded  situations. 

(b)  out-of-plane  stresses  do  not  agree  well,  particularly  near  the  hole. 

(c)  circumferential  strains  agree  well  for  both  loaded  and  unloaded  cases: 
radial  strains  do  not. 

Notice  that,  for  this  case  where  a/b  is  relatively  small,  the  step  changes 
in  o2  and  er  at  the  elasto-plastic  boundaries  are  not  discernible. 

As  noted  above,  as  well  as  the  differing  theoretical  basis  between  the 
analyses,  the  stress-strain  relationships  also  differ  in  the  plastic  region.  It  is 
considered  that,  in  spite  of  these  differences,  the  approximate  theory  given  here 
gives  good  predictions  for  those  parameters  likely  to  be  of  major  interest  -  in-plane 
stresses  and  radial  displacement  (i.e.  circumferential  strain).  For  the  purpose  of 
obtaining  at  least  a  first  estimate,  use  of  the  present  formulations  would  appear  to 
be  worthwhile. 

7.  CONCLUSIONS 

Closed-form  expressions,  which  are  easily  evaluated,  have  been  given  for  the 
quantities  of  interest  associated  with  the  cold  working  of  holes.  Although  these  have 
been  based  on  a  simplified  theory  of  plasticity,  the  results  agree  well  with  those 
from  a  finite  element  analysis  using  the  exact  (incremental)  theory.  The  differences 
which  do  occur  are  either  small  in  magnitude  or  affect  quantities  of  secondary 
interest  only.  The  effect  of  reaming  in  relieving  residual  stresses  is  minimal. 


13.  R.P.  Carey.  Computed  stress  and  strain  distributions  under  interference  fit 
and  after  cold-working.  Dept.  Defence,  Aeronautical  Research  Laboratories, 
Structures  Technical  Memorandum  466,  August  1987. 

14.  R.P.  Carey.  Three  dimensional  computation  of  stress,  strain  and  strain 
energy  density  under  interference-fit  and  after  cold-working  of  holes.  Dept. 
Defence,  Aeronautical  Research  Laboratory,  Structures  Technical 
Memorandum  478,  January  1988. 


(a)  Hole  loaded  by  internal  pressure  (loaded) 


Material  undergoing  plastic 


(b)  As  above  after  removal  of  internal  pressure  (unloaded) 


FIG.  1  IDEALISATION  OF  HOLE  COLD-EXPANSION 


Circumferential 


Ra  dial 


b/a  =  5 
v  =  1/3 


ELASTOPLASTIC  LOADED  ANNULUS 


c/a 


FIG .  4  EFFECT  OF  GEOMETRY  a/b  ON  YIELD  RADIUS  c/a 


FIG.  5  RELATION  BETWEEN  YIELD  AND  REYIELD  RADII  c/a 
AND  p/a  IN  TERMS  OF  GEOMETRY  b/a 
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FIG.  7 


EFFECT  OF  REAMING  ON  RESIDUAL  CIRCUMFERENTIAL 
STRESS 


FIG.  8  STRAINS  IN  ANNULUS  UNDER  LOADED  AND 
UNLOADED  CONDITIONS 
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FIG.  9  CHANGE  IN  STRAIN  RESULTING  FROM  1 0%  REAMING 
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FIG  10(a)  EFFECT  OF  INTERFERENCE  ON  ELASTO-PLASTIC 
RADIUS  (g  =1/3) 


EFFECT  OF  INTERFERENCE  ON  ELASTO-PLASTIC 
RADIUS  (q  =1) 


Stress  (MPa)  Stress  (MPa) 


FIG.  1 1  COMPARISON  OF  IN-PLANE  FINITE  ELEMENT  AND 
ANALYTICAL  STRESSES 


Stress  (MPa)  Stress  (MPa) 


FIG.  1 2  COMPARISON  OF  OUT-OF-PLANE  FINITE  ELEMENT 
AND  ANALYTICAL  STRESSES 


Strain  (x10~3)  Strain  (x  0~3) 


FIG.  13  COMPARISON  OF  FINITE  ELEMENT  AND  ANALYTICAL 
STRAINS 
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